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Abstract 

A particular kind of quintessence is considered, with equation of motion 
Pq/pq = ~~ 1) corresponding to a cosmological term with time-dependence 
A(t) = A(t )(R(t ) I R(t)) p which we examine initially for < P < 3. Energy 
conservation is imposed, as is consistency with big-bang nucleosynthesis, and 
the range of allowed P is thereby much restricted to < P < 0.2. The position 
of the first Doppler peak is computed analytically and the result combined 
with analysis of high-Z supernovae to find how values of Q m and Qa depend 
on P. Some comparison is made to the CMBFAST public code. 
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Our knowledge of the universe has changed dramatically even in the last five years. 
Five years ago the best guess, inspired partially by inflation, for the makeup of the present 
cosmo logical energy density was Q m = 1 and = 0. However, the recent experimental 
data on the cosmic background radiation and the high - Z (Z — red shift) supernovae 
strongly suggest that both guesses were wrong. Firstly Vl rn ~ 0.3 ± 0.1. Second, and more 
surprisingly, Q\ — 0.7 ± 0.2. The value of ^a is especially unexpected for two reasons: it is 
non-zero and it is > 120 orders of magnitude below its "natural" value. 

The fact that the present values of Q m and ^a are of comparable order of magnitude is 
a "cosmic coincidence" if A in the Einstein equation 

R^u — -^QiivR = SttGnT^ + kg^ u 

is constant. Extrapolate the present values of Q m and Vt\ back, say, to redshift Z = 100. 
Suppose for simplicity that the universe is flat Qc — and that the present cosmic parameter 
values are Vl m = 0.300... exactly and f2 A = 0.700... exactly. Then since p m oc R(t)~ 3 (we 
can safely neglect radiation), we find that Q m ~ 0.9999.. and f2 A ^ 0.0000.. at Z — 100. 
At earlier times the ratio Q\/Q m becomes infinitesimal. There is nothing to exclude these 
values but it does introduce a second "flatness" problem because, although we can argue for 
Q m + £7a = 1 from inflation, the comparability of the present values of fl m and Q\ cries out 
for explanation. 

In the present paper we shall consider a specific model of quintessence. In its context we 
shall investigate the position of the first Doppler peak in the Cosmic Microwave Background 
(CMB) analysis using results published by two of us with Rohm earlier [1]. Other works 
on the study of CMB include [2-5]. We shall explain some subtleties of the derivation 
given in [1] that have been raised since its publication mainly because the formula works far 
better than its expected order-of-magnitude accuracy. Data on the CMB have been provided 
recently in [6-13] and especially in [14]. 
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The combination of the information about the first Doppler peak and the complemen- 
tary analysis of the deceleration parameter derived from observations of the high-red-shift 
supernovae [15,16] leads to fairly precise values for the cosmic parameters Q m and We 
shall therefore also investigate the effect of quintessence on the values of these parameters. 

In [1], by studying the geodesies in the post-recombination period a formula was arrived 
at for the position of the first Doppler peak, l\. For example, in the case of a flat universe 
with Vtc = and Qm + Qa = 1 and for a conventional cosmological constant: 



li — 7T 



Rq 



Qm I 



Qa 



V2 „fla 

Rt 



dw 



^/Qmu; 3 + Q\ 



(1) 



If Q c < the formula becomes 



h = 



71 



Rt_ 

\/—Qc 



dw 



(2) 



For the third possibility of a closed universe with Q c > the formula is: 



Tt 



Qc V-Ro 



1/2 



sinh \ Q 



R t 



dw 



i V^m^ 3 + Ha 



(3) 



The use of these formulas gives iso-Zi lines on a Q M — Q A plot in 25 ~ 50% agreement with 
the corresponding results found from computer code. On the insensitivity of Z x to other 
variables, see [17,18]. The derivation of these formulas was given in [1]. Here we add some 
more details. 



The formula for l x was derived from the relation l x = n/ A9 where A9 is the angle subtended 
by the horizon at the end of the recombination transition. Let us consider the Legendre 
integral transform which has as integrand a product of two factors, one is the temperature 
autocorrelation function of the cosmic background radiation and the other factor is a Legen- 
dre polynomial of degree /. The issue is what is the lowest integer / for which the two factors 
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reinforce to create the doppler peak? For small / there is no reinforcement because the hori- 
zon at recombination subtends a small angle about one degree and the CBR fluctuations 
average to zero in the integral of the Legendre transform. At large I the Legendre polynomial 
itself fluctuates with almost equispaced nodes and antinodes. The node-antinode spacing 
over which the Legendre polynomial varies from zero to a local maximum in magnitude is, 
in terms of angle, on average n divided by /. When this angle coincides with the angle 
subtended by the last-scattering horizon, the fluctuations of the two integrand factors are, 
for the first time with increasing /, synchronized and reinforce (constructive interference) 
and the corresponding partial wave coefficient is larger than for slightly smaller or slightly 
larger /. This explains the occurrence of 7r in the equation for the Z x value of the first doppler 
peak written as l\ = n/A9. 

Another detail concerns whether to use the photon or acoustic horizon, where the former 
is v^3 larger than the latter? If we examine the evolution of the recombination transition 
given in [19] the degree of ionization is 99% at 5, 000°K (redshift Z = 1,850) falling to 
1% at 3,000°K (Z = 1,100). One can see qualitatively that during the recombination 
transition the fluctuation can grow. The agreement of the formula for l\, using the photon 
horizon, with experiment shows phenomenologcally that the fluctuation does grow during 
the recombination transition and that is why there is no full factor of v^3, as would arise 
using the acoustic horizon, in its numerator. When we look at the CMBFAST code below, 
we shall find a factor in li of ~ 1.22, intermediate between 1 (optical) and (acoustic). 

To introduce our quintessence model as a time-dependent cosmological term, we start 
from the Einstein equation: 

iV - \ R 9»u = Ht)g^ + BnGT^ = BnGT^ (4) 

where A(t) depends on time as will be specified later and = diag(p, —p, —p, —p). Using 
the Robertson- Walker metric, the '00' component of Eq.(4) is 
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while the 'ii' component is 



'R\ 2 k 8nGp 1 t 

R + W = — + 3 A (5) 



Energy-momentum conservation follows from Eqs.(5,6) because of the Bianchi identity 
D^(R^ - \g^) = D^Ag^ + SttGTV) = D»% v = 0. 

Note that the separation of T^ v into two terms, one involving A(t), as in Eq(4), is not 
meaningful except in a phenomenological sense because of energy conservation. 

In the present cosmic era, denoted by the subscript '0', Eqs.(5,6) become respectively: 

— Po = H S + Jq~ 3 Ao (7) 

-8nGp = -2q H$ + H 2 + A _ Aq (8) 

Kq 

where we have used q = —j^jp and H = 

For the present era, p <C p for cold matter and then Eq.(8) becomes: 

go = t^m - (9) 

where tt M = § jj$ L and fi A = 



Now we can introduce the form of A(t) we shall assume by writing 

A(f) = bR{t)- p (10) 

where b is a constant and the exponent P we shall study for the range < P < 3. This 
motivates the introduction of the new variables 

P ~ 3 

Vt M = Vt M — p^A, ^A = p^A (11) 
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It is unnecessary to redefine Qc because Qm + ^a = ^m + ^a- The case P = 2 was proposed, 
at least for late cosmological epochs, in [20]. 



The equations for the first Doppler peak incorporating the possibility of non-zero P are 
found to be the following modifications of Eqs. (1,2,3). For Vt c — 



1-1=71 



Ik 

R 



V2 „5a 

Rt 



dw 



(12) 



If Qc < the formula becomes 



h = 



7T 



Rt 



-Rn x " 



1/2 



X 



x sin 



dw 



\/Cl M w 3 + Haw p + VL c w 2 , 



(13) 



For the third possibility of a closed universe with Sl^ > the formula is: 



it /R t 



'VLq \Rq 



x sinh \ fl 



1/2 



X 



(14) 



The dependence of l\ on P is illustrated for constant Qm = 0.3 in Fig. 1(a), and for the 
flat case ilc = in Fig. 1(b). For illustration we have varied < P < 3 but as will 
become clear later in the paper (see Fig 3 below) only the much more restricted range 
< P < 0.2 is possible for a fully consistent cosmology when one considers evolution since 
the nucleosynthesis era. 



We have introduced P as a parameter which is real and with < P < 3. For P — > we 
regain the standard cosmological model. But now we must investigate other restrictions 
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already necessary for P before precision cosmological measurements restrict its range even 
further. 



Only for certain P is it possible to extrapolate the cosmology consistently for all < w = 
(Ro/R) < oo. For example, in the flat case Qc = which our universe seems to approximate 
[14], the formula for the expansion rate is 

if (I) =&MW 3 + n A w p (is) 

This is consistent as a cosmology only if the right-hand side has no zero for a real positive 
w — w. The root w is 

If < VLm < 1, consistency requires that P < 3Qm- 



In the more general case of flc ^ the allowed regions of the Qm ~ plot for P = 0, 1, 2 
are displayed in Fig. 2. 

We see from Eq.(16) that if we do violate P < 3il M for the flat case then there is a w > 
where the cosmology undergoes a bounce, with R = and R changing sign. This necessarily 
arises because of the imposition of D^T^ = for energy conservation. For this example it 
occurs in the past for w > 1. The consistency of big bang cosmology back to the time of 
nucleosynthesis implies that our universe has not bounced for any 1 < w < 10 9 . It is also 
possible to construct cosmologies where the bounce occurs in the future! Rewriting Eq.(16) 
in terms of Qa'- 

W = U A -(3-P)J (17) 

If P < 3, then any < will lead to a solution with < w < 1 corresponding to a bounce 
in the future. If P > 3 the condition for a future bounce is < — (^f^)- What this 



means is that for the flat case Qc = with quintessence P > it is possible for the future 
cosmology to be qualitatively similar to a non-quintessence closed universe where R — at 
a finite future time with a subsequent big crunch. 

Another constraint on the cosmological model is provided by nucleosynthesis which re- 
quires that the rate of expansion for very large w does not differ too much from that of the 
standard model. 

The expansion rate for P = coincides for large w with that of the standard model so 
it is sufficient to study the ratio: 

(R/R) 2 P /(R/R) P=0 W =T (3Q M - P)/((3 - P)Q M ) (18) 
«Z5» ( 4 Q R - p)/(( 4 - P)Q R ) (19) 

where the first limit is for matter-domination and the second is for radiation-domination 
(the subscript R refers to radiation). 

The overall change in the expansion rate at the BBN era is therefore 

(R/R) 2 P /(R/R)% =0 W =T (3n M - P)/((3 - P)Q M ) x (40*™ - P)/((4 - P)Q% ans ) (20) 

where the superscript "trans" refers to the transition from radiation domination to matter 
domination. Putting in the values Vl M = 0.3 and Vt l ^ ans = 0.5 leads to P < 0.2 in order that 
the acceleration rate at BBN be within 15% of its value in the standard model, equivalent 
to the contribution to the expansion rate at BBN of one chiral neutrino flavor. 

Thus the constraints of avoiding a bounce (R = 0) in the past, and then requiring consistency 
with BBN leads to < P < 0.2. 

We may now ask how this restricted range of P can effect the extraction of cosmic 
parameters from observations. This demands an accuracy which has fortunately begun 
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to be attained with the Boomerang data [14]. If we choose li = 197 and vary P as 
P = 0,0.05,0.10,0.15,0.20 we find in the enlarged view of Fig 3 that the variation in the 
parameters Qm and Q A can be as large as ±3%. To guide the eye we have added the line 
for deceleration parameter q = —0.5 as suggested by [15,16]. In the next decade, inspired 
by the success of Boomerang (the first paper of true precision cosmology) surely the sum 
(Om + ^a) will be examined at much better than ±1% accuracy, and so variation of the 
exponent of P will provide a useful parametrization of the quintessence alternative to the 
standard cosmological model with constant A. 

Clearly, from the point of view of inflationary cosmology, the precise vanishing of Qc = 
is a crucial test and its confirmation will be facilitated by comparison models such as the 
present one. 

We have also studied the use of the public code CMBFAST [21] and how its normalization 
compares to that in [1]. For example, with P = and fl A = 0.3, h 100 = 0.65 we find using 
CMBFAST that 

n A = 0.5, h = 284 (h = 233 from [1]) 

n A = 0.6, h = 254 (h = 208 from [1]) 

n A = 0.7, h = 222 (h = 182 from [1]) 

n A = 0.8, h = 191 (h = 155 from [1]) 

The CMBFAST Z x values are consistently ~ 1.22 times the li values from [1]. As men- 
tioned earlier, this normalization is intermediate between that for the acoustic horizon (v^) 
and the photon horizon (1). 

Finally, we remark that the quintessence model considered here is in the right direction 
to ameliorate the "age problem" of the universe. Taking the age as 14.5Gy for VL M = 0.3, 
Vtc = and /iioo = 0.65 the age increases monotonically with P. It reaches slightly over 15 
Gy at the highest-allowed value P = 0.2. This behavior is illustrated in Fig 4 which assumes 
Vl M = 0.3 and flatness as P is varied. 
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Figure 1. 

Dependence of l x on P for (a) fixed Vt M = 0.3; (b) fixed VL C = 0. 
Figure 2. 

Regions of the Qm ~ plot where there is a future bounce (small dot lattice), no 
(unshaded) and a past bounce (large dot lattice) for (a) P = 0; (b) P = 1; and (c) P 

Figure 3. 

Enlarged view of £Im — plot to exhibit sensitivity to < P < 0.2. 
Contours are (right to left) P = 0,0.05,0.10,0.15,0.20. 

Figure 4. 

Age of the universe in units 10 w y versus P. 

This figure assumes Vl M = 0.3 and flatness VL C = 0. 
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Figure 1 . (b) 




Figure 2 . (a) 
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Figure 2 . (c) 
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